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The power spectrum of a homogeneous and isotropic stochastic variable, characterised by a finite 
correlation length, does in general not vanish on scales larger than the correlation scale. If the 
variable is a divergence free vector field, we demonstrate that its power spectrum is blue on large 
scales. Accounting for this fact, we compute the gravitational waves induced by an incompressible 
turbulent fluid and by a causal magnetic field present in the early universe. The gravitational wave 
power spectra show common features: they are both blue on large scales, and peak at the correlation 
scale. However, the magnetic field can be treated as a coherent source and it is active for a long 
time. This results in a very effective conversion of magnetic energy in gravitational wave energy at 
horizon crossing. Turbulence instead acts as a source for gravitational waves over a time interval 
much shorter than a Hubble time, and the conversion into gravitational wave energy is much less 
effective. We also derive a strong constraint on the amplitude of a primordial magnetic field when 
the correlation length is much smaller than the horizon. 
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I. INTRODUCTION 

Gravitational waves, once emitted, propagate freely in 
spacetime (affected only by the expansion of the uni- 
verse); they do not interact and are not absorbed. They 
hence provide a direct probe of physical processes that 
took place at the time of their generation, in the very 
early universe. The potential detection of a gravitational 
wave background of primordial origin provides an ob- 
servational signal of sources that have disappeared by 
now. Primordial sources of gravitational waves are for 
example inflation, cosmic topological defects, a pre big 
bang phase of expansion, primordial phase transitions, 
a phase of turbulence in the cosmic fluid and primordial 
magnetic fields. Gravitational waves produced during in- 
flation 0, by topological defects Q or during a pre-big 
bang phase ,3] have been studied extensively. In this pa- 
per we will consider the last two sources which can also 
be associated with a phase transition. 

Gravitational waves are most probably the only rem- 
nant from a period of turbulence in the primordial 
plasma. In order to be maintained, turbulence neces- 
sitates a continuous injection of energy in the absence 
of which it dissipates rapidly leaving no other trace be- 
hind. A possible source for this energy input could come 
from the false vacuum decay during a first order phase 
transition. 

Primordial magnetic fields could be at the origin of 
the ubiquitous magnetic fields observed in galaxies and 
clusters today. Their imprint in the form of gravitational 
waves can confirm their presence in the primordial uni- 
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verse also on scales on which magnetic fields are com- 
pletely dissipated today. 

In the following we study in detail the generation of 
gravitational waves by these two sources and we com- 
pare the results obtained. The sources share common 
properties: they are both related to divergence-free vec- 
tor fields, and they both obey a relativistic equation of 
state so that they evolve like radiation with the expan- 
sion of the universe. They are also both correlated over a 
finite correlation length, if one accounts only for causally 
created magnetic fields and turbulence. However, tur- 
bulence lasts for a short time {e.g. the duration of the 
phase transition) compared to the evolution of the uni- 
verse and it acts incoherently as a source of gravitational 
waves, while the magnetic field is a coherent source ac- 
tive in principle up to the time of matter and radiation 
equality. 

In the next section we demonstrate that the power 
spectra of both sources are blue on scales larger than 
the correlation scale. This implies that turbulence in 
the primordial universe has a Batchelor spectrum on 
large scales. Moreover, it gives rise to a background 
of gravitational waves also on these large scales, which 
has been disregarded so far. In section IIIII we model 
the power spectrum of turbulence and of the turbulent 
anisotropic stress, and we analyse the creation of grav- 
itational waves. We propose a new method to perform 
this analysis, avoiding an error in the gravitational wave 
dispersion relation which is found in previous works. We 
find that the gravitational wave spectrum peaks at a fre- 
quency corresponding to the size of the largest turbulent 
eddies. In Section lTVl we compute the gravitational waves 
induced by a causal magnetic field correlated on a scale 
which may be much smaller than the horizon, and we 
constrain the magnetic fields using the nucleosynthesis 
bound. Finally, in Section[V] we compare the two results, 
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concluding that a radiation-like, long acting, coherent 
source like the magnetic field is more efficient in gener- 
ating gravitational waves on large scales than a rapidly 
dissipating, incoherent source such as turbulence. We 
believe that this result is valid beyond the two examples 
discussed here. 

We consider a Friedmann universe with flat spatial sec- 
tions, 

ds 2 = a 2 (r])[-dr] 2 + S^d^dx 3 } . (1) 

The variables x, r, L etc. denote comoving distances and 
k is the comoving wave vector. The density parameter 
is always scaled to today, Clxiw) = Px(v)/ Pc(vo), where 
the index o indicates the present time. For relativistic 
species we have therefore Qx(v) = ^x( f lo)/ a4 ( 1 l) : > we 
normalise a (770) = 1 and sometimes denote the present 
value of a density parameter simply by flx(vo) = ^x; 
likewise, p c = Pc(Vo)- 7~L = f denotes the conformal 
Hubble parameter. The radiation energy density today 
is taken to be Q ra d(f7o) = ^rad = 4.2 x 10~ 5 0. We use 
the expression 

for the scale factor which is strictly true only in a flat 
matter-radiation universe. But since we are mainly con- 
cerned about the physics during the radiation dominated 
phase of the universe, neglecting the cosmological con- 
stant A is irrelevant for our purpose. Values of confor- 
mal time T] <C 2(\/2 — l)V^rad/-ffo = ?7cq denote the ra- 
diation dominated phase, while during the matter domi- 
nated phase 7] 3> 77eq. 

II. THE POWER SPECTRUM OF 
DIVERGENCE FREE VECTOR FIELDS AT 
LARGE SCALES 

We want to discuss gravitational waves generated from 
the energy momentum tensor of divergence-free vector 
fields. In this section we derive the generic form of the 
spectrum for such vector fields. 

We consider a stochastic vector field v(x), statistically 
homogeneous and isotropic, hence («i(x)) = 0. The gen- 
eral form of its two point correlation tensor is (see for 
example @) 1 

M x ) = My My + x)> = 

= T,(r)(5ij — XiXj) +T(r)xiXj , (2) 

where x is the unit vector in direction x, and r is the 
length of x. Here E(r) and T(r) denote the correlation 
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under parity (see 0). We neglect parity violating terms in this 
work. 



functions in directions perpendicular and parallel to x. 
If the stochastic vector field is divergence-free (for ex- 
ample in the case of the turbulent velocity field of an 
incompressible fluid, or in the case of a magnetic field), 
the correlation functions S(r) and T(r) are related. It is 
easy to verify that 3 g'j x ^ = yields 

^ = ;(=(r)-r(r)). (3) 
dr r 

Homogeneity implies that the Fourier transform of v for 
different wave vectors k ^ q is uncorrelated, so that the 
power spectrum is of the form 

(ViCk)t$(q)) = (2^) 3 <5(k - q)P y (k) , (4) 

where Py(k) = J dPse* - "^ (x) . (5) 

Isotropy and V • v = require the following simple form 
for the power spectrum Py (k) 

P ij (k) = P(k)(5 ij -k i k j ). (6) 
The trace of Eq. (0 together with Eq. J5J lead to 

P(k) = J d 3 :re* k ' x (V) + ir(r)) . (7) 

In the expanding (non inflationary) universe, any causal 
random process is characterised by a finite correlation 
length L, which is at most the size of the horizon at the 
epoch at which the process is occurring, L < 77*. The 
correlation function bij (x) is generally expected to decay 
at large distances, but in the case of the expanding uni- 
verse, causality even requires that it vanishes identically 
on scales larger than the horizon rj* (if t!j(x) is a clas- 
sical field, see 0). Therefore, in the case of interest to 
us, the correlation function 6y(x) of a causal process is 
always a function with compact support. This implies 
that its Fourier transform must be an analytic function 
in the variable k. Hence, the power spectrum © must 
be analytic in k. In order for P(k){5ij — k{kj) to be an 
analytic function of the components ki also at k = 0, we 
have to impose P{k — > 0) oc k n , with n an even positive 
integer, n > 2. From Eqs. J^J and © it follows that 
P(k) is given by 

P{k) = \ J d 3 x^b u (r)=2n [°° drr 2 ^lb u (r). 

(8) 

For small k we obtain up to order k 2 

P(k -> 0) = 2tt / drr 2 b vl (r) 
Jo 

poo 

- -k 2 drr 4 6 l2 (r)+0(fc 4 ). (9) 
Jo 

Analyticity of Pij (k) requires the first term on the right 
hand side of Eq. @ to vanish. This result can also be 
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obtained directly from Eq. @ using the fact that T(r) 
vanishes for sufficiently large distances r. Indeed, using 
(0 one finds 



drr 2 bn(r) =r 3 r(r)|~ =0, 



(10) 



since T(r) has compact support. One is thus left with 



P(k^0)=(^-J^ drr 4 T(r)^ k 2 +C(fc 4 ). (11) 

These considerations apply both to the power spec- 
trum of a magnetic field, if it has been generated in the 
early universe by a causal process Q, and to the power 
spectrum of turbulent motions in the primordial fluid, 
which can be generated for example at the end of a first 
order phase transition (if one assumes that the fluid is 
incompressible) U El El El 



III. GRAVITATIONAL WAVES FROM 
TURBULENCE 

A period of turbulence in the primordial universe can 
leave an observable trace today in the form of a stochas- 
tic background of gravitational waves. The most natural 
source of turbulence in the early universe is a first order 
phase transition (see EH f° r another potential source of 
turbulence), where bubbles of the broken phase nucle- 
ate and expand into the symmetric phase. Therefore we 
concentrate on this case. We parameterise th e sp ecifics 
of the phase transition following 0, 0, 0, Il2| . The 
temperature and the value of the scale factor at which it 
takes place are denoted by T», a*; the comoving size of 
the largest bubbles at the moment when they collide is 
L. This length corresponds to the largest scale on which 
turbulence develops (the size of the largest eddies). The 
ratio of the false vacuum energy density to the thermal 
energy density of the universe at the moment of the tran- 
sition is a = Pvac/ P*'i and the velocity of the bubble walls 
Vb(a). Finally, the fraction of vacuum energy density 
which is converted into kinetic energy, turbulence, of the 
fluid is denoted by k = pidn/Pvac- 

To generate a considerable amount of gravitational 
waves, one needs to move big masses rapidly. This can 
be achieved when a first order phase transition proceeds 
through detonation of the broken phase bubbles. In this 
case the bubble walls propagate faster than the speed of 
sound 1/V3 < Vb(a) < 1, and there is a large concentra- 
tion of energy at the bubble walls. The phase transition 
is completed in a (conformal) time interval of the or- 
der of A77* ~ L/vb(a) which is much smaller than one 
Hubble time Ar]*H* <C 1. Typical values for Ar/^H* 
for the electroweak phase transition are 0.001 — 0.01 
(see 0, 0, 0, an d references therein) . 

A first order phase transition causes injection of ki- 
netic energy in the primordial fluid at a characteristic 
scale L corresponding to the size of the largest bubbles 



at the moment when they collide. Therefore, the pri- 
mordial plasma is stirred on the scale L <C Ti^ ■ Since 
the Reynolds number is very high in the early universe, 
turbulent motions are generated EH- We assume that 
the generated turbulence is homogeneous, isotropic, fully 
developed and stationary. Then, since L is the principal 
turbulence scale, on smaller scales (larger wave numbers) 
k > L^ 1 an energy cascade is established which pro- 
ceeds from large to small scales. Eventually one expects 
the formation of a Kolmogorov spectrum in the 'inertial 
range' L -1 <fc< A - , where A denotes the dissipation 
scale In the next section we derive an approximate 
form of the turbulent power spectrum which is valid at all 
wavenumbers (including large scales, small wavenumbers 
k<L- x ). 



A. Turbulence power spectrum 

The variable w^(x) represents the fluid velocity at the 
position x ] and we denote its power spectrum by P v (k). 
Following 10] we define the total kinetic energy per unit 
enthalpy of the fluid as 

_^L-=i< Ui (x) Vi (x)) = ifc ii (0) ) (12) 
p* + V* 2 2 

where p* and p* denote the energy density and pres- 
sure at the time of the turbulence. This non-relativistic 
approximation is strictly valid only for (v 2 ) <C 1: the 
Kolmogorov theory of turbulence has been formulated 
and verified only in this regime. However, we shall use 
it also in the relativistic regime, imposing the upper 
bound (v 2 ) < 1/3 (EH) where it still gives the cor- 
rect order of magnitude. The details of the Kolmogorov 
spectral shape are probably not valid in the relativistic 
regime, but we follow here earlier literature on the sub- 
ject, see 0,E3,E1E]. As usual, we also define the energy 
dissipated per unit time and unit enthalpy by 0, llH Il7j| 



£ = — - 



dr) 2 



(13) 



The turbulent velocity power spectrum can be evalu- 
ated from Eq. (JSJ) , if the shape of the two point correla- 
tion function bu(r) is known. If the turbulence is homo- 
geneous and isotropic, in the inertial range A <C r <C L 
the correlation function is given by the Kolmogorov two 
thirds law (er) 2 / 3 . In the dissipation range r -C A, the 
fluid motions are regular and one can Taylor expand the 
velocity field as a function of r. This entails that the 
correlation function is proportional to r 2 Moreover, 
the symmetric tensor field fry(x) has compact support 
and its divergence vanishes. On this basis, we make the 
following Ansatz for the transverse correlation function 




for < r < A 

for A < r < L 
otherwise , 



(14) 
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with Eo an as yet unspecified constant. From the 
divergence-free condition @ we now determine the ra- 
dial correlation function 



r(r) 



3 





.2/3 



for < r < A 

for A < r < L 
otherwise . 



Eq. J2J also requires that T be continuous, so that E does 
not contain a singularity (r may, however, have a kink 
which leads to a discontinuity in E). The continuity of 
r(r) fixes E in terms of A, (v 2 ) and L: T(r > L) = 
implies 



2/3 



(15) 



(where we have neglected a term of order (A/L) 2 ). The 
radial correlation function is non negative, continuous 
and identically zero on scales larger than the correla- 
tion scale. The transverse correlation function starts 
positive, passes through zero inside the inertial range at 
r = (3/4) 3/2 L, and has a jump of -(l/9)(v 2 ) at r = L in 
order to satisfy the condition of compact support. This 
jump corresponds to a kink in T(r) for r = L. The in- 
tegral of bu(r) over space vanishes, so that the power 
spectrum satisfies P v (k = 0) = 0. Fourier transforming 
the correlation function given above, one finds 



P v (k) ~ 2tt(v 2 ) x 



51.2 



TO.", 



E°k 



for 



§ x/3r(§)L- 2 /3A;-ii/3 for L- 1 « k « A" 1 (16) 



81 





otherwise . 



In the inertial range we recover the Kolmogorov spec- 
trum, but the power is non-zero also for small wave- 
numbers k < L~ x (a well known result in the theory 
of turbulence [Til Il7| , but which has not been taken into 
account in cosmological applications [ToL ITlj 2 so far). 

The discontinuity of the spectrum at the peak k = L^ 1 
is not real. It is a consequence of our approximations 
which are valid only for kL <C 1 and kL 1 respec- 
tively. A numerical calculation of the spectrum would 
smoothly turn from the k 2 behavior for kL <C 1 to the 
f or kL 1. The discontinuity indicates, how- 
ever, that the peak amplitude obtained by extrapolating 



2 In Refs. IH H3, |2<| |3(1 the spectrum at large scales 
is given, but there MHD turbulence, not pure fluid turbulence 
is considered. Furthermore, the fact that for a divergence free 
vector field causality implies that the s pect rum is blue, P v oc k 2 
for Lk -C 1 is only acknowledged in Ref . |24| . If we a ctually follow 
the renormalization group analysis outlined in 1 3 1 1 - we find that 
at large scales, kL <C 1, the driving force has a k 4 spectrum 
and the stable fixed point has the critical exponents z = 2 and 
X = —3/2. This can be inferred simply from causality and from 
the fact that the field is divergence free. 



the asymptotic behavior right up to k = 1/L is not very 
reliable. Eq. l|12f> gives the spectrum normalisation 



MO) 



d 3 k 

(27T) 5 



Pv(k). 



(17) 



Integrating Eq. (|l(jfl . we obtain instead 0.76(i; 2 ). This 
overestimation of the total energy is also due the exten- 
sion of our approximations up to the peak k = L^ 1 . It 
indicates that the peak amplitude might be roughly a 
factor of 2 too high. In the following, we will neglect 
these discrepancies which are within the accuracy of our 
approximations. 

We define the ratio of the total kinetic energy of tur- 
bulence to the radiation energy density (c/. Eq. Ijl2|l ): 



fir fa*) _ 1 



^rad fa* 



&ia.d(V*) JO 
jOkin _ 2(t> 2 ) 
P* 3 



dk dflr{k, 77*) 
k dlog(fc) 
d 3 k 



6tt 3 



P v (k). (18) 



The above equation defines the turbulent energy spec- 
trum per logarithmic unit of frequency 



dClT(k, 77*) 
dlog(fc) 



ftradfa*) T^k 3 P v (k) . (19) 



Usually the turbulent energy spectrum is defined as 
E(k) — k 2 P v (k)/(2Tr 2 ) p|, in our notations we have 



(20) 



In the Kolmogorov inertial range the energy spectrum 
has the form E{k) = Ce 2 ^k~ 5 ^ 0, following H3 we 
set the value C = 1. Comparing E(k) and the power 
spectrum (|16|) we find the relation 



2{eL) 



2/3 



(21) 



On large scales k < L^ 1 the above considerations on 
the properties of null divergence and compact support of 
the correlation function imply the formation of a Batch- 
elor spectrum for the turbulent energy. The definition of 
the energy spectrum <|19|) and Eq. 1111) give 



dfln 



dlog(fc) 



^radfa*) 7T Ik 
97T 



(22) 



where I = J °° dr r 4 T(r) is Loitsyansky's integral |j| 0, 
[l7| . Since the correlation function has compact support, 
the Navier-Stokes equation requires that I be constant 
in time. Therefore, the turbulence is persistent on large 
scales and its rate of decay increases with wavenumber. 
Loitsyansky's integral is a measure of the angular mo- 
mentum of the turbulent fluid. Its constancy is a man- 
ifestation of the conservation of angular momentum . 
A simple estimate gives I ~ (v 2 )L 5 . Exploiting the con- 

d{v 2 ) ^ jv 3 ) 
dr) ~ L ' 



stancy of I and the energy decay equation 
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Kolmogorov has estimated the decay of the kinetic energy 
in freely evolving turbulence (v 2 )(rj) oc ?y _10 / 7 and the 
growth of the correlation length L(rf) oc rj 2 / 7 [M. fill ll7T|. 
In some cases these decay laws are a reasonable fit to 
the experimental data, nevertheless nothing tells us that 
they are verified for our oversimplified model of turbu- 
lence in the primordial fluid. In the following, we will 
simply neglect the overall decay of turbulence as far as 
gravitational wave production is concerned. We shall as- 
sume that turbulence appears at some time 77,,, remains 
active for a short period r and simply disappears after- 
wards. We shall not model its decay in any detail. 



B. Eddy turnover and the time of duration of 
turbulence 



The energy cascade involves a hierarchy of vortices of 
different sizes. To each of these structures on a given 
scale one can associate a characteristic velocity, the eddy 
velocity. Following Refs. \WL lllj. we define the charac- 
teristic eddy velocity at the largest scale L by 



2 ^rad(??*) 



dk d£lT(k, 
k dlog(fc) 



3 nrfa.) 

2 firad(?7*) 



(23) 



where the second equality holds since most of the turbu- 
lent energy resides in the large wavelengths k ~ L^ 1 . The 
velocity of the largest eddies, v^, is another parameter 
that measures the energy stored in turbulence. For the 
maximal eddy velocity, v\ — 1/3 we have q T ^* \ = §■ 
In this regime, our non-relativistic approximation under- 
estimates the kinetic energy of turbulence somewhat. 

The velocity of the eddies on a length scale A < £ < L 
is simply 



VL ~~L 



1/3 



(24) 



One further defines the eddy turnover time ti = £/vi and 
the turnover frequency u>£ — 1/t£. 

This picture is valid only on scales X < £ < L. The 
principal scale of turbulence, L, defines the extent of the 
region over which velocities are correlated, and it is in- 
terpreted as the maximal size of turbulent eddies. On 
scales larger than L there is no correlated fluid motion 
and therefore no eddies. Similarly on scales smaller than 
A the Reynolds number is of order unity and the fluid is 
not in a turbulent regime, therefore no eddies are formed, 
and we do not associate a turnover frequency on scales 
larger than L or smaller than A. However, the energy 
distribution is not zero for £ < A and £ > L since the 
eddies of a given size £ contribute to the power spectrum 



P v (k) across the full range of wavenumbers 3 |l7|. 

The eddy turnover time on a given scale defines the 
time interval over which the eddies of that size break 
down into smaller ones, transferring energy from the orig- 
inal scale to a smaller one. This is the principle of the 
direct cascade. If the stirring of the fluid is continuous, 
fully developed and stationary turbulence is established 
in a time interval of the order of the turnover time on the 
largest scale tl- Once the stirring stops, the turbulence 
decays away in the same time interval. If the stirring 
lasts for less than t£ , the cascade of energy develops but 
the turbulence is not stationary. In Ref. [10( it is argued 
that, from the point of view of gravitational wave pro- 
duction, this case can be modelled like fully developed 
stationary turbulence lasting for a time tl (instead of 
for the stirring time). The assumption which goes into 
this model is that eddies of every size £ < L are gener- 
ated and dissipated in the time interval tx (so that they 
make several oscillations before transferring energy to a 
smaller scale). Therefore, we follow ^(| and define the 
total time of duration of turbulence as r = max(Tx, Ary*). 
We take here Arj^ because the external stirring lasts until 
the phase transition is completed. 

For a first order phase transition this always reduces 
to r = t l 10] . From Eqs. and J2D we have 



1 L 2 / 3 
TL ~——. (25) 



y^e 1 ' 3 

Since pkin = p* = p*/3, Eq. l(T3*)l gives 



3 net 



(26) 



Let us assume that the stirring lasts for much longer than 
the turnover time of the largest eddies. Then Ary* should 
be substituted in the above equation for e, so that 



tl 



J_ /4\ 1/3 L 2 ' 3 ^) 1 ' 3 
y/2 W (««) 1/3 



(27) 



With this expression and using that L ~ Vb(a)Ar]^, the 
condition A^* 3> tl can be rewritten in terms of the 
phase transition parameters, 



A77* > T L 



Vb{a) 3k 

^ <<: VT 



Vb(a) and n(a) are given in Ref. [lTJ: 



v b (a) 
n(a) 




a" + -a 



(28) 

(29) 
(30) 



3 In the previous treatment (c/. Eq. 1161 ). we have neglected power 
on scales smaller than A. The dissipation scale contains very 
little energy, and is too small to be relevant in the context of a 
cosmological treatment. 
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Substituting these expressions in (|28() it turns out that 
the condition A77* 3> T£ is never verified for a < 1. 
Therefore, even in the case of a strongly first order phase 
transition, the turbulence persists for a duration of the 
order of the turnover time of the largest eddies. In the 
following analysis we set r ~ tl- During this time we 
neglect the decay of the turbulent kinetic energy and 
the subsequent time dependence of the power spectrum 
P v (k) . We assume that the turbulence is stationary dur- 
ing the time interval tj,, and instantaneously dissipated 
just afterwards. 



C. Anisotropic stress power spectrum 

We now evaluate the anisotropic stress generated by 
the turbulence. 

The energy momentum tensor of the turbulent fluid 
has an anisotropic stress component that induces tensor 
perturbations hij of the metric. The propagation equa- 
tion for tensor perturbations in a radiation dominated 
universe is 

2 • 

hij(k, 77) + -hij(]s.,rf) + k 2 hij{k,rf) = 87rGa 2 (r?)n y (k, 77) , 

V . , (31) 
where Ily (k, 77) is the anisotropic stress. Tensor pertur- 
bations are induced on each mode k. If the mode is super- 
horizon at the moment of generation, the last term on the 
left hand side of the above equation is much smaller than 
the damping term -hij and the tensor perturbation is not 
yet oscillating. Once the perturbation enters the horizon, 
the oscillatory term can no longer be neglected and the 
perturbation becomes a gravitational wave in the proper 
sense. 

If the source of the anisotropic stress has a blue power 
spectrum on large scales k < L^ 1 , then the anisotropic 
stress power spectrum is white noise (and not zero) for 
every k < L^ 1 as we now show. At the present time 
the modes Hq < k < L^ 1 are sub-horizon. On these 
scales there is a background of gravitational waves com- 
ing from the primordial source, characterised by a white 
noise spectrum. In the case of turbulence after a phase 
transition, this signal has not been calculated in previous 
analysis. 

The energy momentum tensor of the fluid is 



Tpv = (p + p) VfiVv +pg fiU - 



We extract the spatial, transverse and traceless part 
which is relevant for tensor perturbations. The tur- 
bulence is active on sub-horizon scales k I/77* and 
for a time much shorter than one Hubble time. For 
this short period of time, the expansion of the uni- 
verse can be neglected, and we may assume a Minkowski 
background. In Eq. (|32() only the second order term 
Iij(x, 77*) = 4/3 p^,Vi(x, r}*)vj(x, 77*) generates tensor 
perturbations |1C| . We have substituted = p*/3 since 
the phase transition takes place during the radiation 



dominated era. In k-space, the tensor anisotropic stress 
is given by %(k) = (P a P 3m - l/2P ij P lm )T lm (k) where 
Pij = 5ij — kikj is the transverse traceless projector, and 
Tij(k,r)*) = fp* / ^^«i(q,7?*)'Wj(k - q,??*), up to an 
irrelevant trace which vanishes in the projection. The 
anisotropic stress power spectrum (IL^ (k, 77* )n*„- (q ,77*)) 
has already been calculated in Ref. [lOj: using Wick's 
theorem to reduce the four point spectral function of the 
velocity, one obtains 

(n ij (k ) r7*)n*,(q,r7*))= (^A S(k - q)Il v (k, t,.) , 



n„(fc,j7*) 



(33) 



d 3 pP4p,77*)P.(|k-p|,77*)(l+ 7 2 )(l+/3 2 ), 



, (34) 

with 7 = k • p, = k • k — p. In principle, we are not al- 
lowed to use Wick's theorem, the turbulent velocity field 
being not at all Gaussian. We invoke it here for simplic- 
ity and to obtain an expression for the four-point func- 
tion. We hope that even if not correct in detail, this gives 
the right order of magnitude for the four point function. 
Reducing the four point spectral function to the power 
spectrum by means of Wick's theorem is sometimes used 
as a closure method also in more detailed studies of tur- 
bulence [13 ■ 

The anisotropic stress power spectrum is explicitly cal- 
culated at the time of the phase transition 77, . The time 
evolution equation for the turbulent power spectrum is 
highly non linear and involves the three point function [j| . 
As already mentioned in Section^ in our approximate 
approach to the evaluation of gravitational waves gen- 
erated from the turbulence, we ignore the detailed time 
dependence of the anisotropic stress, since it is not very 
relevant for the resulting gravitational waves which are 
given by a time integral of the source. This however does 
not mean that we are treating the turbulence as a coher- 
ent source, as will become clear in the next section. 

The integral in l|34l) cannot be done exactly, but good 
approximations have already been derived in Ref. |lfj for 
small scales k > L~ l , and in Refs. 0,H(j for large scales 
k < Lr x . With Eqs. (QUl and JTSJl, we find 4 



( 32 ) \l v {k) 



27 



8 \^rad(f7*) 



2 r3 

for k < L 1 

L -2/3 fc -ll/3 

for L^ 1 < k < A -1 
. otherwise . 



(35) 



4 The analytical approximations used in the above mentioned ref- 
erences are not entirely correct, because of the bounds that have 
been used in the angular integration. 7 = k p is not varying from 
— 1 to 1 but its bounds depend on k and p. We have checked 
that using the correct bounds does not significantly modify the 
approximated final result. 
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As anticipated, the tensor source power spectrum is white 
noise on large scales. However, because of the convolu- 
tion, the largest contribution to n„(A;) on these scales 
does not come from the velocity power spectrum on the 
same scales, but from the value of the velocity power 
spectrum at the maximum, k ~ L^ 1 . The discontinuity 
of the spectrum at k = L^ 1 is not physical, it is due to 
the approximations we used in evaluating integral • 



D. Generation of gravitational waves 

We now determine the spectrum of gravitational waves 
which arc induced by the turbulent motion of the primor- 
dial fluid. The turbulent velocity is a random variable, 
but within Richardson's model of the energy cascade it is 
possible to associate a frequency to the fluid motions at 
each scale, since the fluid motions are vortical. This 'os- 
cillatory' behaviour is peculiar to each realisation of tur- 
bulence, and is lost in the statistical average. The veloc- 
ity power spectrum P v (k) is not oscillating in time. First 
of all, the Fourier transform decomposes the velocity field 
into a hierarchy of waves, not of eddies. Moreover, a sta- 
tistical model cannot reproduce the detailed dynamics of 
individual vortex events. However, the 'oscillatory' be- 
haviour of the source is relevant for the generation of the 
gravitational wave background. 

In order to account for these oscillations which are 
present in turbulent motion on scales smaller than L, we 
model the random velocity of the turbulent fluid heuris- 
tically as 



wi(k) 



for k < L~ 



v i {k)e iu » r > for L" 1 < k < A" 1 



(36) 



where we do not associate a turnover frequency to scales 
larger than the correlation scale. Here lj^ is the frequency 
of an eddy with size £ < L, 

vt/e = K/£ 1/3 )r 2 / 3 

(v L /L l l 3 )k 2 ' 3 « k, k= l/t . (37) 



The energy momentum tensor of the source in k space 
involves a convolution (cf. section UlI C|l : 

TV (k, = J ^« ^(cO^fk-qJe^e*-!-- 1" . 

The convolution entails that the energy momentum ten- 
sor as a random variable is oscillating also on large scales 
k < L^ 1 , contrary to the random velocity field. At these 
scales the amplitude is dominated by the contributions 
from q — l/L (cf. Eq. (|35fl ). We therefore expect the en- 
ergy momentum tensor to oscillate with the smallest fre- 
quency This can be approximatively justified by the 
following reasoning. Since we concentrate on the case of 
turbulence generated by a first order phase transition, the 
source is active for a time tl- Therefore, the main con- 
tribution to the convolution integral comes from frequen- 
cies such that bJqTL ~ 1, W| k _ q |Ti ~ 1. If k < L^ 1 and 



q < L^ 1 the frequencies are both zero, but for k < L^ 1 
and q > L~ x one can approximate wik— q| ~ w <? so that 
the frequency is 2u> q . The same argument for smaller 
scales k > L^ 1 gives the frequency 2ujk- We can therefore 
model the time dependence of the random anisotropic 
stress tensor by 



IIy-(k,7?) = 



__ J nij(k)e 2iULr i for k < L^ 1 



n y (k)e 2,w *'' for L- 1 <k < X- 1 



(39) 



With this source the gravitational wave equation be- 
comes the equation of a forced harmonic oscillator with 
damping, 

2 • , 
hij ■ (k, i]) + - h l3 ■ (k, i]) + k 2 hij (k, ti) = (40) 

= 8^Ga 2 (77)n 4 ,(k) e 2 ^e(77 - Vin )Q( mn - r,) . 

Here the source frequency is generically denoted 2u>, and 
rji n and ?7fin are the initial and final times of action of 
the source, 7yg n — r/i n = tl- The Heaviside functions 
are inserted to limit the source to the interval r/j n < r\ < 
77fi n . The source is active for a time much shorter than 
one Hubble time: T£ <C ?/*. Therefore we can neglect 
the expansion of the universe during the relatively short 
period during which the source is active, 



hij (k, T]) + k 2 hij (k, 77) = 



(41) 



87rG^n y -(k)e 2iu " ) e(ry - 7^)6(7^ - 77) . 



The solutions to the homogeneous equation without 
damping are simply exp(±ifc?7), and with the Green's 
function method one finds the inhomogeneous solution 

for 77 < 77 fln , 



hijQz, 77 < 77 fin ) 



AnGalRijik) 



Jkr) 



i(2w— k)r} _ i(2a)-fc)?; in 



k - 2uj 



-ik7] & ^ 

k + 2us 



(42) 



Note that the singularity at k — 2uo of the first term 
is removable. For r/in < i] < ?7 nn the time dependence 
of gravitational waves with wave number k is given by 
the superposition of the source frequency, 2u>, and of the 
source wavenumber k. However, once the source ceases 
to be active, the gravitational wave propagates freely at 
the speed of light. In order to find a solution for r\ 3> 77 nn , 
we need to match Eq. I|42|l at 77 nn with a generic solution 
of the homogeneous propagation equation 

2 • , 
hj (k, 77) + -hij (k, 77) + k 2 h^ (k, 77) = . (43) 

Now the expansion of the universe can no longer be ne- 
glected, and the homogeneous solutions in this case are 
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exp(±ikr/)/r). The matching procedure leads to 
47rGa^*II ii (k) 



e ikr] e i(2ui—k)r]t in _ p i(2u—k)rn n 



7] k — 2uj 

g— ikri e i(2ui+k)r] Cia _ gi(2w+fc)r;;, 



V 



k + 2ui 



(44) 



(where we have substituted rygn by 77* in the amplitude). 
For 77 ^> 77* the gravitational wave propagates freely with 
the correct dispersion relation: its frequency is just ±k, 
and the time dependence r? -1 accounts for the expansion 
of the universe. The frequency 2w of the source only 
enters in the amplitude. 

We want to evaluate the gravitational wave energy den- 
sity today, normalised to the critical density p c . In real 
space the gravitational wave energy density parameter 

is m 



G = ^ = 

Pc 



!6TrGp c a 2 



(45) 



Here (• • •) denotes ensemble average and, at the same 
time, time averaging over several periods of oscillation. 
The factor 1/a 2 comes, since the overdot denotes the 
derivative w.r.t. conformal time. 

We define the gravitational wave power spectrum by 



(MM)/C-(q,77)> = S(k- q)|h| 2 (M) 



(46) 



Fourier transforming the above relation for £l G one de- 
fines the gravitational wave energy density per logarith- 
mic unit of frequency: 



dk dQ G (k) 
k dlog(fc) 



with 



k 3 \h\ 



dn G {k) _ 

dlog(fc) ~ 4(27r) 6 Gp c a 2 



(47) 



(48) 



This is the quantity we are interested in. To proceed we 
need to find an expression for the spectrum 146(1 at late 
times (today) and on sub- horizon scales krj 1. From 
solution (|44|l with kr\ 3> 1 one has, up to an irrelevant 
overall phase 



MM) 



47rGa 2 77*IIy (k) 



^4 e ife (>J-')in) _ £ e — ik(»J-»Jta) 



(49) 



where A 



and B 



i(2u, + k)T L _ 



Given the 



k-2ui fc+ 2o; 

anisotropic stress power spectrum l|33|) . the gravitational 
wave power spectrum takes the form 



\h\\k,ri) = ^GfaH (^ 2 n t ,(fc,^) |A|2 + |i?|2 . 

V (50) 

We now use Eqs. l(1%|l and (f5T7f> to evaluate the gravita- 
tional wave energy density spectrum. The main spectral 
dependence comes from the anisotropic stress spectrum 
((35(1 . which we write in the form 



JL 1-3 
13 

for < x < 1 

x -2/3 (51) 

for 1 < x < j 
, otherwise , 



rn„(fe) ~ — 7T 4 



27 4 / n T {v*) 



8 vO ra a(j?*) 



where we have defined the variable x = kL. Moreover, 
the value of u> in the amplitudes A and B depends on 
k. From Eq. ((39(1 . we have u> = lol for scales k < L~ l , 
and lo — LOk <x A: 2 / 3 for scales fc > L _1 . Using the ex- 
pression for the scale factor in the radiation dominated 
era a — iJov^rad??, we finally obtain the energy density 
spectrum, 



dQ G (k,r)) 



dlog(fc) 64tt 



9 iWfW" 



^rad(?7*)/ a 4 (r?) 



rad 



13 



(x-2v L ) 2 



"L -1 



(:E-2t; L iE 2 / 3 ) 2 



(x+2 t ) i ) 2 



2 "i -1 



(x+2?;!,:!; 2 / 3 ) 2 



for < x < 1 



2.-2/3 for 1 < x < j 



otherwise . 



(52) 



Once the source has decayed, the gravity wave energy spectrum today is simply given by the above expression 
density evolves like radiation, as it should. The energy evaluated at 770 with 0(770) = 1. We obtain a different 
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: i — i — i — i — r_r 

10" 6 - 





10" 8 r 



10 - 9 r 

0.1 1 10 

kL 



FIG. 1: The gravitational wave energy density today, nor- FIG. 2: Like Fig. □ but with vl = 1/2, hence ± 
malised to (Tts,L) 2 fi ra d, as a function of x = kL is shown. 



»t(i.) 



1/6. 



10' 



The black, oscillat- 



We have chosen vf = | ijrfa*) 
ing line shows the exact amplitude given in Eq. 15211 . The 
red, straight line shows the approximated spectrum given in 
Eq. I53L normalised to the same quantity (7i*L) 2 f2 ra( j. The 
factor 2/13 appearing in Eq. (1521 has been neglected in the 
plot since it generates an unphysical discontinuity. 



spectral dependence whether the largest eddy velocity 
vl is below or above 1/2. For vl < 1/2, the characteris- 
tic frequency 2u>l is less than the scale L _1 at which the 
anisotropic stress power spectrum peaks (|51ll . The grav- 
itational wave energy spectrum still peaks at but it 
changes slope at k — 2lol- An approximate solution for 
the energy spectrum in the case vl < 1/2 is 



dn G (k, m ) 9 2 
dlog(fc) ~ 32^ { * ' 



^radfa) 



rad 



X < 



x 3 /v L for < x < 2vl 

4 x for 2vl < x < 1 

4a;~ 8/3 for 1 < K j 

, otherwise . 



(53) 



The value of the energy spectrum normalised to the pa- 
rameters (H*L) 2 f2 ra d is shown in Fig.^for vl = 0.01 and 
in Fig. H for v L = 1/2. We remind that = \v\ 

(cf. Eq. m- 

If v L > 1/2, then 2uj k > 2u L > L~ x . The velocity of 
the largest eddies is bounded by the speed of sound vl < 
1/V3 0, 0. Also in this case, the energy spectrum 
peaks at the frequency L —1 , but because of the particular 
form of the amplitude it changes slope at k = (2vl) 3 /L. 



An approximate solution for vl > 1/2 is 



rf» G (Mo) _ J_ (r/ n 2 / 
dlog(fc) " 32^ [n * ' \ 
x 3 jv\ for < x < 1 



' »r fa) 

v ^rad(?7*) 



ad 



X < 



X' 



■~ 2 ' jv\ for 1 < x < (2v L ) 3 
4s" 8 / 3 for (2v L ) 3 < x < j 
K otherwise . 



(54) 



This result is plotted in Fig. for the maximal value 

v L = l/\/3, ft^h = I" Eqs - ® and ® expressed 
in terms of the parameters which describe the phase tran- 
sition (see Refs. H3, E3), 

can be found in Appendix B. 

The maximal gravitational wave signal is obtained for 
the highest possible value of the eddy velocity, vl- Given 
that L ~ Vh{a)lS.ri^, we obtain a peak frequency of 



8 x 10" 3 



1 



1 



T, 



v b {a) Ar]*H* 100 GeV 



mHz (55) 



(we normalise it to the electroweak phase transition tem- 
perature). Moreover, the phase transition has to be 
strongly first order, so that the bubbles expand at the 
speed of light. If Vb(a) ~ 1 and Ar]*H* ~ 10 -2 as ar- 
gued before, the spectrum peaks at a frequency of about 
mHz. For the maximal eddy velocity of vl = we 
obtain the amplitude at the peak frequency 



dQ G {k,r] ) 



dlog(fc) 



-J- x HT 4 fi rad ~ 10- 10 . (56) 

Z47T 



The total energy density in gravitational waves for 1/2 < 



10 




FIG. 3: The figure shows the gravitational wave energy 
density today, normalised to (7i»L) 2 f2 ra( j, as a function of 
x — kL, for vl = lA/3, hence ~ T ^*\ = I . The red, 
straight line is the approximation given in Eq. 154H . 

vl < l/v3 is given by the integral of Eq. I|54|) . 

Jo k dlog(fc) 
*sf-(-Y£%Vw (57) 
For the maximal value v L = 1/V3 we have therefore 

— . . [ J ^rad 

144tt \W 
~ 0.01f h 2 (a) O rad . (58) 

We find that the turbulent phase arising after a strongly 
first order phase transition produces a background of 
gravitational waves which is in principle detectable by 
the space interferometer LISA. However, the gravity wave 
intensity obtained here depends strongly on the peak am- 
plitude of the turbulent spectrum, which is rather uncer- 
tain. Therefore, also the amount of gravitational waves 
generated is probably not very accurate. 

Our result is in agreement with Ref. ^H, with a few 
differences. Since we account for the natural dispersion 
relation of gravitational waves u) = k, we obtain that 
the spectrum peaks at the frequency 1 /L, and not at the 
turnover frequency lu l as found in 9, 10, 111. The two 
differ by a factor v L . In our case the characteristic size of 
the eddies determines the frequency of the induced gravi- 
tational waves, while the eddy velocity enters only in the 
overall amplitude (see Appendix A). We therefore also 
have a different slope of the spectrum for k > L~ x than 
Refs. [9T fToL ITiT | . Moreover, we find that the gravitational 



wave spectrum is not reduced only to the 'Kolmogorov' 
part, but it continues to wavelengths larger than 1/L 
as shown in Figs ^ to For a strongly first order phase 
transition, the amplitude of the signal at mHz frequencies 
derived here is of the same order of magnitude as found in 
|llj |. whithin the range of our analytical approximations. 

IV. GRAVITATIONAL WAVES FROM 
MAGNETIC FIELDS 

A. The magnetic field power spectrum 

We now analyse the production of gravitational waves 
by a primordial magnetic field. In order to make a di- 
rect comparison with the turbulence case, we consider 
here only causally created magnetic fields. Typically, 
like turbulence, magnetic fields can be generated dur- 
ing a phase transition (note however that in the case of 
the electroweak phase transition, magnetic seed fields can 
form even if the phase transition is second order Q). 

The magnetic field power spectrum on large scales is 
determined along the same lines as the turbulence power 
spectrum in Section [H] The field is divergence- free and 
it must be uncorrelated above a certain scale, therefore 
the slope has to be k n with n > 2 an even integer. 
This slope continues up to the correlation scale. If the 
magnetic field is simply constant on scales smaller than 
the correlation scale the power spectrum goes to zero 
steeply on these scales (in 24] for example it is mod- 
eled as exp(-(/c//c c ) 4 )). This is normally assumed to be 
the initial spectrum. However as the system evolves, the 
interactions of the magnetic field with the cosmic fluid 
modify the spectrum at small scales, turning it into a 
power law (for example, into the Iroshnikov Kraichnan 
spectrum, Ea(k) = C(ey AY^ 2 k~ 3 ^ 2 with va the Alfven 
speed, see El El 13 To stay general, we make the 
ansatz 

(B t (k)B*{q)) = (27r) 3 <5(k - q)(% - fc. t ^)P B (fc) , (59) 

( C B k 2 for 0<k<L- 1 

P B (k) = I C B L a - 2 k a for L- 1 <k< X- 1 (60) 
I otherwise , 

with a < —3, so that the integrated energy is dominated 
by the contribution from the correlation scale L. L and 
A are not necessarily the same scales as for turbulence 
but are the magnetic correlation and dissipation scales. 
Cb is an as yet undetermined constant. 

Once the magnetic field is generated, it is not quickly 
dissipated like turbulence. Because of the high conduc- 
tivity of the cosmic plasma there is negligible Ohmic dis- 
sipation and the magnetic flux is conserved. During the 
early stages of evolution of the universe both, the kinetic 
and magnetic Reynolds numbers are very high, and mag- 
netohydrodynamical turbulence is generated. The mag- 
netic field decay has been analysed in this case both ana- 
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lytically 26, 27j,|2J,|29j and numerically |2J,|25j, and dif- 
ferent scaling laws are obtained depending on the treat- 
ment of the problem. In the following we keep the scaling 
exponent as a free parameter, since it will not be very rel- 
evant for the final result. However, all papers agree on 
the fact that the power spectrum is persistent on large 
scales, i.e. that Cb is constant in time. 

After electron-positron annihilation the neutrinos are 
already decoupled and the photon mean free path in- 
creases suddenly . Then the system is no longer turbu- 
lent but dissipative due to radiation viscosity, and the 
magnetic field evolution changes [22I I23L |25| . 

The main point for our analysis is that in the process of 
generation of gravitational waves, one needs in principle 
to account for the evolution of the magnetic spectrum. 
This can influence the final result, since the magnetic 
field is acting as a source of gravitational waves for a 
cosmologically long time. The situation is different from 
the case of turbulence. This evolution has not been taken 
into account in |Ts| , where a fixed correlation length given 
by the size of the horizon at the time of generation has 
been assumed. 

As for the slopes of the magnetic spectrum (|6L)fl . we 
make a general Ansatz also for its time evolution. The 
evolution is mainly determined by three processes: the 
expansion of the universe, the growth of the correlation 
length and the decay of the energy density. To eliminate 
the time dependence by redshifting we work with the 
comoving magnetic energy density 

(B i ) = (B 3 (f,))a\r } ). (61) 

However, as we shall see below, (B 2 ) is still time depen- 
dent due to dissipation of energy. 

As in the case of the turbulence, we normalise the spec- 
trum l|5U)l in terms of the total magnetic energy density 
divided by the radiation energy density: 

n B (v) = 1 (b 2 ) 

^va,d(v) Pra,d(Vo) 87T 

1 1 f°° dk 3 

« « J T k ' PB{k > 
dk dils 



which gives 



C B 



Pt3 
1 

^rad Jo k dlogk 



(2tt) 3 5(a+3) 2 5 
8tt a-2 [ ' 



(62) 



(63) 



Here we have defined the magnetic field energy density 
parameter Qb which corresponds to the magnetic field 
energy density scaled to today via Eq. iJBTJ. Qb is itself 
a function of time due to the dissipation of energy. 

Let us now study the time dependence of L and {B 2 }. 
We model the comoving correlation length by 



m = m - 



with < 7 < 1 so that L{rj) never overcomes the horizon 
at any instant r\. L* denotes the correlation length at 
the epoch 77* of creation of the magnetic field. The time 
dependence of the dissipation scale A does not play a 
significant role in our analysis, so we do not consider it. 
For our purposes, it is enough to have in mind that A is 
also a function growing with time. In order to find the 
time evolution of the magnetic energy density, we make 
an analogy: we assume that from the MHD equation it is 
possible to define a conserved quantity analogous to the 
Loitsyansky's integral. We impose simply 



(B 2 )L 5 — constant in time, 



(65) 



with this the large scale part of the spectrum is persistent 
(see Eq. (|63|) above) . This is equivalent to the constancy 
of Cb , which has been observed in simulations [3 l2. l: 
With this the comoving magnetic energy scales as 



(B 2 )( V ) = {B 2 )M[^ 



57 



(66) 



If we were to adopt the same decay laws as in turbulence, 
one would find the Kolmogorov decay laws L(r]) oc rj 2 ^ ', 
(B 2 ){ri) oc rf^l 1 . Equivalent scaling laws are obtained 
from the arguments of self similarity used in |24ll26j . and 
of selective decay used in [25, 28] . The main difference is 
that most of these works consider a white noise spectrum 
for the magnetic field at large scales, resulting in a slower 
decay of the energy and a more substantial growth of the 
correlation length. They find Saffman's law for fluid dy- 
namics. However, as demonstrated in Section^] causally 
generated, divergence-free vector fields in the early uni- 
verse necessarily lead to the formation of a Batchelor 
spectrum. This is valid for both, turbulence and mag- 
netic fields. 



B. Power spectrum of the anisotropic stress 



From the (comoving) magnetic field stress tensor 



T (B) = _j_ 



l -B 2 



— BjB n 



(67) 



(64) 



one defines again the anisotropic stress as n^k) = 
(PuPjm — l/2PijPim)Ti m (k). One then calculates the 
comoving anisotropic stress power spectrum in the usual 
way, see e.g. pj| , using Wick's theorem to reduce 4-point 
to 2-point correlators under the assumption of Gaussian- 
ity (which again is most probably not strictly valid but 
can be expected to give the right order of magnitude). 
Setting 

(Ikj (k, 77)n?. (q, 77)) = n B (fc, T05(k - q) (68) 
one obtains, exactly like in Eq. (|34JI ■ 

n B (M)= / dVM^)^(|k-pU)(i + 7 2 )(i + /3 2 )- 

(69) 
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A good approximation to this convolution integral is 



where we can identify 



f A\L-\r{) 
for < k < L 1 



II B (fc,77)-^H 



A\ L a - J (r)) k c 



2 ; a-7 / 

for < fc < A -1 
otherwise , 



(70) 



with A\ 



7(2a+ 4 3) aild ^ 2 ^ V 5(2q+3) ' 



ct-2 



The time evolutions of the magnetic energy density 
and of the correlation scale are deterministic. The same 
holds for the evolution of the anisotropic stress power 
spectrum. We therefore may treat magnetic fields as a 
coherent source, i.e. a source with deterministic time 
dependence, 

(n iJ (k, % )I^(q,» ft )) = 
/(^i)/(^)(n y -(k,^ in )n* J (q,, ?in ))^k - q) , (71) 

{f(rj) denotes here an arbitrary function). Using 
Ii B (k,77) = / 2 (r/)n B (k,?7 in ), we find 



(II tf (k,T&)IIJ J .(q J r&)) = y/lL B (k,r h Wa B (k,r l 2)S(k-q) . 

(72) 



C. Generation of gravitational waves 

In the case of a coherent source, one can write the evo- 
lution equation (|31|l directly for the gravitational wave 
power spectrum. To see this, we first write it for the 
stochastic comoving anisotropic stress of the magnetic 
field 



2 • 

Tj 



k 2 hu = 8ttG^ 
a 2 



(73) 



The solution to this equation with initial condition 

ha (k, v ) = SrrG J n dQ G(v, C, k) , (74) 

where Q(rj, £, k) is the corresponding Green function. 
The gravitational wave energy power spectrum becomes 



d(- 



d V Jr,, fl2 (C) 



x ^ f d ^a4f {UlJ (k ' mh (q ' 0) • (?5) 
For a coherent source, using l|72ll . we find 

(htjMhiifav)) - N 2 (M)^(k-q) = (76) 

d 



«5(k-q) 



8ttG 



d( ToTT^ yWC,k) 



2 (C) 



7tj. a (0 



(77) 



Since hij(k,r]) and h(k,rj) share the same Green's func- 
tion, h(k,rj) satisfies equation (|73|l . with source term 
87rGy / ns(fc, i])/a 2 . In the radiation era this is 



-h + k 2 h 



8ttG 



-v5b(M) 



77 Hln TaA rf 
With the power spectrum 1(75)1. the solution is 



(78) 



h(k,rj) 



4ttG i 

T 



e -ikr) gifef 



ife?7 p -ifeC 



»7 C »7 < 
for 0< fc < L _1 (0 

A2jL (a-7)/2 (C) fcQ /2 

for L- X (C) < * < ^(0 



(79) 



I 



otherwise . 



For a given mode k, the time integration has to be divided 
in different ranges: if the mode satisfies initially k < L~\ 
then there exists a time ?7z,(fc) at which the correlation 
length overcomes fc -1 , and we have k > L _1 (r/) for 77 > 
rjL^k). From this time on, one has to consider the small 
scale part of the source spectrum. Moreover, the modes 
that are super-horizon at beginning fc < 777 1 have to 
be considered separately, since in this case the integrand 
does not oscillate until horizon crossing. However, since 
the magnetic correlation length L(rj) is increasing with 
time, the above integrals are in any case dominated by 
their values at the lower bound 77*. We split the solution 
in three fc regions, and we finally obtain 



h(k,r)) 



4ttG i 

— Gb 



f^rad fc 
A T -7/2 2 

for < fc < f]^ 1 



A X L 



-7/2 



for 77- 1 < fc < Ll 1 

a 2 L: 7/2 (kL*r/ 2 ^ 

for L~ x < k < A" 1 



(80) 



'/ 



I 



otherwise . 



This solution has to be matched with the solution of the 
homogeneous equation at the time the source ceases to 
be active: that is, for every mode fc, the time ri\(k) at 
which the dissipation scale has grown bigger than fc -1 . 
However, the above solution does not depend on ri\(k), 
being dominated by the lower bound 77* , and has the same 
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time dependence as a free propagating wave. Therefore, 
the above solution is simply valid also for 77 > r]\(k). 

We want to evaluate the gravitational wave energy den- 
sity per logarithmic unit of frequency given in Eq. 1)48 [I. 
For this we have to average \h\ 2 over several periods. 
Substituting the expression of the scale factor a(rj) = 
i/oV^rad^ we find 

iii.( M) _as£c. 2 



z £ rad 



a 2 (77) 



for < k < ri^, 1 

> 1 1 

for i]~ L < k < L~ L 



(81) 

AlL-T{kLX{ikY fOT L ^ < k < K 1 








otherwise . 



0.00010.001 0.01 0.1 

L,k 



10 100 



Here we have used that h ~ kh on sub-horizon scales. 
The gravitational wave energy density is now given by 
(|48p. We substitute the constant Cg using expression 
pty. and since HbL 5 is time independent, we can insert 
its value at the initial time, rJs(?7*)L^. These manipula- 
tions lead to (see fig. 0J 



dS~l G {k,r]o) 
dlog(fc) 



rad 



A? (^) 2 (£*fc) 3 for < k < v- 1 
A\ {^)\Uk) for 77- 1 < fc < L" 1 



(82) 



A\ (^Y(L,kr+ 1 for L- 1 < k < A- 1 
otherwise . 
The total energy density can be approximated by 
dk dflc(k, 770) 



Or 



75tt 2 



a-2 



dlog(fc) 
2 

A? 



(83) 



1 



2 ^(77*) 



Or. 



ail 



where we have neglected a term (L^/i]*) 3 , given that 
L* < r/*. If the magnetic field at small scales devel- 
ops an Iroshnikov Kraichnan spectrum corresponding to 
a = —7/2, the above quantity becomes 



77 V7* 



fi|(?7*) 



rad 



(84) 



D. Limits on the magnetic field intensity 

In paper [Tsj . we found strong constraints on the am- 
plitude of a primordial magnetic field by applying the nu- 
cleosynthesis bound on the gravitational waves generated 
by the magnetic field. The limits apply to magnetic fields 
created both during inflation and during a phase transi- 
tion, and are formulated in terms of the magnetic field 




10 100 



FIG. 4: The gravitational wave energy density today, nor- 
malised to f2s/S! ra d, as a function of x = fcL» is shown for 
L,TL* = 0.01 in the top panel and for L*TL, — 1 in the bottom 
panel. On scales k < I/77* the spectrum behaves as fc 3 , then 
it flattens into linear growth up to the correlation scale L~ x . 
On smaller scales it decays like k a+1 . In this plot a — —7/2 
and 7 = 2/7 are chosen. 



intensity Ba present today on the scale A = 0.1 Mpc. 
This is the relevant scale for the fields observed in clus- 
ters of galaxies. Since the magnetic field generation takes 
place long before 77 = 0.1 Mpc ~ 10 13 sec, the mode A is 
super-horizon at the moment of creation. 

If the magnetic field is generated at inflation, it is char- 
acterised by a simple power law spectrum k n with n > 
—3, up to an upper cutoff which corresponds to the time 
at which inflation ends, k < rj^ 1 , where 77* ~ 8 x 10 -9 sec 
is a typical value. The generation of magnetic fields dur- 
ing a phase transition instead is a causal process, leading 
mainly to fields on scales smaller than the size of the 
horizon at the moment the field is created. In the case 
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of the electroweak phase transition at 100 GeV, this cor- 
responds to 77* ~ 10 5 sec. 

In Ref. 0| w e have neglected sub-horizon fields, which 
cannot propagate to larger scales unless they are helical 
and subject to an inverse cascade. Motivated by infla- 
tion, we have imposed a cutoff scale 77" 1 in the magnetic 
power spectrum corresponding to the horizon size at the 
phase transition. We have argued, that this leads to con- 
servative bounds for the magnetic fields. Here we show 
that taking into account also the sub-horizon contribu- 
tion strengthens the limits somewhat. 

As we have seen, gravitational wave production takes 
place also inside the horizon. Moreover, under realistic 
circumstances, most of the magnetic energy is stored on 
these scales. Therefore, we now recalculate the bounds 
obtained in [lsf accounting also for sub-horizon modes. 
The results we find below reduce to those of Ref. [18| in 
the limiting case L» ~ 77*. As for the previous sections, 
our results reside on the assumption that the large scale 
part of the magnetic spectrum is persistent. 

The nucleosynthesis limit on any additional radiation- 
like form of energy comes from the fact that this latter 
may not significantly change the expansion law of the 
universe during nucleosynthesis. The maximum allowed 
additional energy density redshifted to today is about 
;$ 10~ 5 |21|. We impose this bound on the gravi- 
tational wave energy density generated by the magnetic 
field. To be specific, we now fix the value a = —7/2, but 
this choice does not affect the final result in a relevant 
way. From Eq. I|84|l . we have therefore 



0|(?7*) 



< 



77 

24^2 



at I 



(85) 



The upper bound on the total magnetic energy density is 
increased for L* < 77*. We want to formulate the bound 
in terms of the magnetic field amplitude on the scale 
A = 0.1 Mpc. At this aim, we perform a volume average 
of the field in a region of size A 3 . We convolve the field 
with a Gaussian window function 



B A (x) = J ^B(x + y)exp 



y 

2A 2 



(86) 



and we define the magnetic field energy density smoothed 
on a scale A, B\ = (Ba(x) • Ba(x)). We have therefore 



2 1 

n .10 



dkk 2 Ps(fc) exp 



A 2 fc 2 



(87) 



Since A » L„ the above integral is dominated by the 
large scale part of the spectrum fcL* < 1, and inserting 
Eq. lt)3|) we obtain finally 



r2 120 nr a t n \ ( L * 



were again we have taken flsL 5 at the initial time. We 
can now use Eq. (|85|l to derive an upper bound on the 



magnetic field intensity B\: 
B A 



lO" 6 Gauss 



< 



ty) \l) 



5/2 



(89) 



Comparing with Eq. (33) of Ref. [HI, we see that the 
bound is strengthened by a factor (^p) ■ The reason for 
that is twofold. First, the limit on Qb is increased by 
a factor r/*/L* when reducing L» for fixed 77*. On the 
other hand, the magnetic field is peaked at the smaller 
scale L„: this appears in the value of B\ through the 
factor (L*/A) 5 = 5 (77* /A) 5 . 

In the case of the electroweak phase transition at 
100 GeV, for which 77* ~ 10 5 sec, we find 



B A < 10-^ I hi 

V* 



Gauss . 



(90) 



If the phase transition is first order and proceeds via bub- 
ble nucleation, we expect the correlation length to be 
given by the size of the largest bubbles which has been 
estimated as L* ~ 0.017/,, leading to the gravity wave 
spectrum shown in the top panel of Fig. However, if 
the phase transition is second order, the magnetic field 
inside one horizon volume is created smoothly and has 
enough time to align itself. In this case we therefore ex- 
pect a correlation length of the order of L* — leading 
to the gravity wave spectrum shown in the bottom panel 
of Fig. H 



V. COMPARISON OF THE TWO CASES 

We have calculated gravity wave production from two 
sources which are both represented by divergence free 
vector fields. One of them (the magnetic field) is per- 
sistent over many Hubble times while turbulence decays 
within a fraction of one Hubble time. We would natu- 
rally expect the former to produce more gravity waves 
than the latter. However, this is not properly the case. 
Here we compare the efficiency of the two sources, tur- 
bulent and magnetic, in detail. Let us first compare the 
total gravitational wave energy produced, as a function 
of the total energy available in the sources. For this, we 
consider Eqs. I|57|) and l|84|). where the first one applies 
for the largest eddy velocity 1/2 < vl < l/v3 and the 
second one for the value of the magnetic field spectrum 
a = —7/2. Denoting the quantities concerning turbu- 
lence by a superscript T and those concerning magnetic 
fields by B, we have 



0.02 



( »r(r|.) \ 



(91) 



The gravitational wave energy density generated by mag- 
netic fields depends quadratically on the total energy 
of the magnetic source, while for turbulence the rela- 
tion is linear. This is due to the presence of the factor 
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v l 2 = f%^77 m amplitude of the gravitational 
wave turbulent spectrum (|54|l . If the turbulent and mag- 
netic correlation scales as well as the total energy avail- 
able in the sources are comparable, the efficiency in gen- 
erating gravitational waves is similar. For vl = l/v3, 
we find from Eq. (|91|l Hq/Hq ~ 0.1, roughly the same 
order of magnitude within our accuracy. We find no sig- 
nificant domination of the magnetic field induced total 
gravity wave energy. 

We now compare the gravitational wave energy den- 
sity spectra, in order to analyse the efficiency of the 
sources as a function of wavenumber k. We want to relate 
Eqs. I|54|) and l]82p. For this we set the correlation lenghts 
to the same value Lj = Lf = L». Therefore we take 
~ 0.01, as in the case of turbulence generated af- 
ter a phase transition (c/. section ITTlj) . Moreover, we fur- 
ther impose that the initial available energy densities are 
comparable, O r (77* )/17 rad (77* ) ~ n B (jj*)/O rad . Applying 
the reasonable values vl = 1/V3, a — —7/2, 7 = 2/7, 
the ratio of Eqs. (|54|l and l|82l) becomes (x = fcl/*) 



x 2 for L*7Y» < x < 1 
0.1 { ^-x 1 ' 2 for 1< x < {2v L f (92) 



rflog(fc) 
/ rfng(fc,r) ) \ 
V dlog(fc) J 
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, otherwise . 



for (2v L ) 3 < x < L*/A* 



From the above expression we can conclude that, on 
super-horizon wavelenghts x < L/H*, k < ?j~ , the 
energy density of gravitational waves from turbulence 



is smaller than the one from magnetic fields by a fac- 
tor 0.1(L*H*) 2 ~ 10~ 5 . For sub-horizon modes, the 
ratio is increasing and reaches the maximum value at 
kmax — (2wl) 3 /-/j*, where the two amplitudes are com- 

P arable ^^(kf / ^dt^ - °- L Since the spectrum 
peaks at k = ~ /s m ax, also the total energy densi- 

ties are comparable. 

We hence find that on super-horizon scales magnetic 
fields are more efficient than turbulence in generating 
gravitational waves. This is due to the difference of the 
time intervals over which the sources act, and to the 
characteristic way in which a radiation-like source con- 
verts energy into gravitational waves. In the integrals of 
Eq. (O, the Green's function oscillates with frequency 
k and decays as a power law. Therefore, super-horizon 
modes, krj < 1 are more efficiently converted into gravita- 
tional waves than sub- horizon modes. Since the magnetic 
field is active for many Hubble times (approximately from 
its generation up the epoch of matter-radiation equality) , 
the magnetic super-horizon modes generate gravity waves 
until they cross the horizon (subsequent generation can 
be neglected). Turbulence, on the other hand, is active 
as a source of gravitational waves only for a time inter- 
val much shorter than one Hubble time tl <C Ti.* 1 . Most 
turbulent modes that start super-horizon do not have the 
time to convert their energy in gravitational waves. This 
leads to the suppression factor (L*/^*) 2 = (L*H*) 2 for 
the gravity wave spectrum from turbulence. 

The expression for the turbulence to magnetic field ra- 
tio in full generality (with — L T — Lf) is 



rff2g(fc,?) ) 
dlog(fc) 

dJ2g(/c,jj ) 
dlog(fc) 



800tt 3 V a + 3 



1 \n Tad (r U )J 

v 2 ( n B {, h ) \ 2 
V n rad J 



' ^ (^) 2 {L*H*) 2 for < x < LM* 
jp x 2 for L*H* < x < 1 
^ x -(a+S) f or 1 < x < (2 VL f 

a .-(a+ii/3) for ( 2vl ) 3 < x< L./A, 



(93) 



otherwise . 



VI. CONCLUSIONS 

We have analysed the generation of a gravitational 
wave background by a period of turbulence in the pri- 
mordial universe, arising after a strongly first order phase 
transition, and by the presence of a primordial magnetic 
field. Since the early universe is an ionised plasma, treat- 
ing pure Kolmogorov turbulence on its own is proba- 
bly not realistic. However, given the previous literature 
on this subject, the analysis in Section IIIII is essential 
EH LLjl LLJ ■ We have present the more realistic case of 



MHD turbulence in Section HVI and a comparison of the 
two cases in Section Ivl 

The turbulent velocity field and the (causally cre- 
ated) magnetic field are necessarily uncorrelated on scales 
larger than a given correlation scale, which is at best 
equal to the cosmological horizon. We have demonstrated 
that this fact, together with the property of vanishing 
divergence shared by the two fields, implies that their 
power spectra are blue on scales larger than the correla- 
tion scale. This was already known in the case of causal 
magnetic fields [32|. but is new for turbulence. A phase of 
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turbulence in the primordial universe is inevitably char- 
acterised by a Batchelor energy spectrum. A Saffman 
spectrum is excluded, since it relies on long range corre- 
lations that can be generated by the pressure fields [Tt| . 
On small scales, we have assumed Kolmogorov two thirds 
law for the turbulent velocity correlation, and we have de- 
rived the turbulent spectrum from direct integration of 
the correlation function. We have used the asymptotic 
behavior of the turbulence power spectrum right up to 
the peak region. Therefore, the amplitude of our result 
is probably not very reliable. We plan to improve this 
with a numerical analysis |33j . 

Since turbulence is an incoherent source, the analysis 
of the gravitational waves generated by it is not straight- 
forward. In Section IlII Dl we propose a procedure where 
the time dependence of the turbulent flow is not lost in 
the ensemble averaging. This time dependence is impor- 
tant for the correct determination of the induced gravi- 
tational waves. Our method to analyse the formation of 
gravitational waves is based on an heuristic model of the 
oscillations of the random velocity field. This model is 
inspired by Richardson's energy cascade through eddies 
of different sizes. 

Our work corrects an error in the dispersion relation 
found in previous papers HE IH 111 and provides the 
correct dispersion relation for the free propagation of the 
gravitational waves , af ter the turbulent source ceases to 
be active. In Ref. [10| the sub-horizon power spectrum 
of the gravitational waves, at the time when the turbu- 
lent source ceases to be active, is given in Eq. (44) as 
a function of k. The component (hij (k, r]*)h*j(c[, 77*)) oc 
<5 3 (k — q) of the spectrum has momentum k. However, 
then they set the gravitational wave frequency at time 77* 
to ujki the oscillation frequency of the turbulent source, 

Uk = K/i 1/3 )fc 2/3 ■ (94) 

This means that the gravitational wave dispersion rela- 
tion at time 77* is the same as the dispersion relation of 
the turbulent source. However, for 77 > 77* the gravita- 
tional wave propagates freely, and the dispersion relation 
should simply be co — k. For 7/ > 77* the gravitational 
wave power spectrum changes only by redshifting with 
the expansion of the universe. Therefore, in this pro- 
cedure it is not clear how one can recover the correct 
dispersion relation for gravity waves at 77 > 77* . 

The origin of this mistake is the idea that the grav- 
ity wave frequency is determined by the frequency and 
not by the wave number of the source. This is certainly 
true for spatially small sources with infinite lifetime |34j . 
However, in Appendix A we show that for a gravitational 
wave source which is localised in time but not in space, 
the frequency of the emitted gravitational wave is de- 
termined by the wavenumber of the source (the generic 
relation between the frequency of the gravity wave and 
the frequency and the wave number of its source are de- 
rived in Ref. |35|). This is particularly important for 
turbulence, since the turbulent eddies do not spin at the 
speed of light. As a consequence, we find that the energy 



power spectrum of the induced gravitational waves peaks 
at a frequency corresponding to the size of the largest 
turbulent eddies (the correlation scale), and not to their 
frequency. The spectrum continues to both larger and 
smaller wavelenghts with different slopes. 

Usually, the nucleating bubbles of true vacuum also 
generate anisotropic stresses in the fluid and therefore 
lead to the production of gravity waves. The spectrum 
of these gravity waves from the bubbles also peaks at 
frequency L~ l . Its amplitude has been estimated for 
example in Ref. 0- 

In the case of magnetic fields, we obtain a blue power 
spectrum on large scales, leaving as free parameters both 
the slope of the power spectrum on small scales, and the 
growth of the correlation length with time. Assuming the 
existence of the invariant quantity (B 2 )L 5 analougous 
to Loitsyansky's integral, we obtain that the comoving 
anisotropic stress depends on time only via the magnetic 
correlation length. We have solved the evolution equation 
for gravitational waves sourced by magnetic anisotropic 
stresses accounting for the evolution of the magnetic cor- 
relation length. This has allowed us to derive new limits 
for the amplitude of the magnetic field as a function of 
the correlation scale, extending the results of Ref. 0] to 
correlation scales which are smaller than the horizon. 

We conclude our analysis with a comparison of the 
efficiency of turbulent and magnetic sources. When im- 
posing that the total energies available in the sources and 
their correlation lenghts be comparable, we find that the 
total gravitational wave energy densities are of the same 
order, but the intensity of gravitational waves generated 
by magnetic fields on super-horizon scales is enhanced 
with respect to that of gravitational waves from turbu- 
lence by a factor (L*H*) 2 . We can conclude that a co- 
herent, relativistic, long acting source as the magnetic 
field is more efficient in generating gravitational waves 
on super-horizon scales, than a source which acts only 
over a brief period of time like turbulence. 
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APPENDIX A 

In this appendix we want to clarify that, for a stochas- 
tic gravitational wave source which is homogeneously dis- 
tributed in space and which has a time duration com- 
parable to its characteristic oscillation period, the fre- 
quency of the generated gravity wave is connected to the 
wavevector k (hence to the typical size) of the source and 
not to its frequency. The calculation is already performed 
for the particular case of the turbulence in section IlII Dl 
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but it is repeated here in more generality for clarification. 
We neglect cosmological expansion since it is of no rele- 
vance for our argument and consider simply a Minkowski 
background. The wave equation 

(d 2 - A)hij (x, t) = 87rGHy (x, t) , (Al) 

can be solved by the retarded potential 

hu (x, t) = SvrG / dV %( x/ ^-| x - x D _ (A2) 
J l x - x l 

For a long lived source confined in some small spatial 
region, one may consider the usual wave-zone approx- 
imation. In the wave zone, that is at distances much 
larger than the dimension of the source, the field emit- 
ted by the source looks like an outgoing spherical wave 
with wavcnumbcr k = uj, where uj is the source frequency 

n«(x,t) = n ij (x)e iuf ||. 

However, we are treating the case of a stochastic, sta- 
tistically homogeneous cosmological source. The source 
is not stationary, but characterized by a well defined time 
interval of action. Therefore, in order to solve the above 
equation we first Fourier transform it in space 

{d 2 + k 2 )h i3 (k, t) = BrrGIIy (k, t) . (A3) 

The retarded Green's function of the above operator (the 
Fourier transform of the Green function used in eq. (|A2|) ) 
is simply 

Gfl(fc,t) = 0(t)^M j (A4) 

where Q(t) is the Heaviside function. If the source starts 
emitting at a time t m we have 

&«(k,t) = 8ttG [ rft , Sm(fc( *~*' )) II i3 -(k,0 . (A5) 

Jti n % 

At the time fg n , the gravitational wave is given by 

MMU = **gJ^ d t' sin{k{tfi » ~ ^ nijOc.o 

^•(k,f fln ) = 8^G^" f,n *'cos(fc(t fin -i'))ny(k,t')- 

After the source stops emitting, the gravitational 
wave background obeys the free wave equation (df + 
k 2 )hij(k,t) — 0, with solution 

, /, x 4 sin(fc(i — tfin)) „ /, / 

% (k, t > *&,) = Ay V V fc hD " + Bij cos(fc(t - t fin )) . 

(A6) 

Matching the free wave and the sourced wave at time t nn 
we obtain 

A (J (k) = 8vrG y *'cos(fc(t fin -t'))n i3 (k,i') 
By(k) = SvrG dt' Sln(fc(tfi fc n ~ t0) n i3 - (k, t') (A7) 



Regardless of the particular time dependence of the 
source ITj(k, t'), the free gravitational field l|A6|l oscil- 
lates with the frequency uj — k corresponding to the 
wavenumber of the source. If the source is itself os- 
cillating, but it obeys a non-trivial dispersion relation 
a; (A;) ^ k, then the source frequency intervenes in the 
amplitude of the free propagating wave. Only if the du- 
ration tfi n — Un is much larger than the period l/u>, the 
above integrals can be approximated by S (k— u>) and the 
wave inherits the frequency of the source |35j . 

The same result can also be obtained from Eq. I|A2(1 . 
using that |x — x'| has to lie in the interval [t— tg n , t — t- m ]. 

This has been misrepresented in previous literature 
(see e.g. 0, EI)- It is not very important for the re- 
sults obtained as long as vl — 0(1), but it is an interest- 
ing conceptual point: A stochastic spatially homogeneous 
source for gravity waves which has a lifetime comparable 
to its characteristic oscillation period (as it is the case 
for turbulence) imprints its wave number and not its fre- 
quency in the gravity waves it produces. On the other 
hand, a stationary and isolated source always imprints 
its frequency and not its wave number. 



APPENDIX B 



In this appendix we rewrite the spectra of the gravi- 
tational waves generated by turbulence in terms of the 
phase transition parameters as they are convention ally 
defined. We conform to the notations of Refs. [Til Il2| 
which are commonly used in the literature. We denote 
us = vl the velocity of the largest eddies, so that the 
total kinetic energy of the turbulence is n T ^"\ — |uo- 
The duration of the phase transition is /3 _1 = a* A??*, 
and is the physical Hubble parameter at the phase 
transition. The frequency of the gravitational waves is 
denoted / = k, and /<? = uj^/a^ is the physical turnover 
frequency of the largest eddies scaled to today. The peak 
frequency given in Eq. Ij55(l corresponds to kh — fs/us- 
We call /max = 1/A the maximal frequency correspond- 
ing to the dissipation scale. Finally, the gravitational 
wave energy density per logarithmic unit of frequency is 
denoted Qcw(f)- We now express our results in terms 
of fs, u s , #*, [3 and v b . 

Instead of considering the full expression (|52|l . we find 
it more useful to rewrite the approximated formulas l|53|) 
and In the case us < 1/2, Eq. we find 
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In the case us > 1/2, Eq. I|54|l. we find 
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